Combined first-principles and EXAFS study of structural instability in BaZrOg 
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Phonon spectrum of cubic barium zirconate is calculated from first principles using the density 
functional theory. Unstable phonon mode with the R25 symmetry in the phonon spectrum indicates 
an instability of the cubic structure with respect to rotations of the oxygen octahedra. It is shown 
that the ground-state structure of the crystal is M/mcm. In order to find the manifestations of the 
predicted instability, EXAFS measurements at the Ba Lm-edge are used to study the local structure 
of BaZrOa at 300 K. An enhanced value of the Debye-Waller factor for the Ba-0 atomic pair 
[al ~ 0.015 A^) revealed in the experiment is associated with the predicted structural instability. 
The average amplitude of the thermal octahedra rotation estimated from the measured af value is 
~4 degrees at 300 K. The closeness of the calculated energies of different distorted phases resulting 
from the condensation of the R25 mode suggests a possible structural glass formation in BaZrOs 
when lowering temperature, which explains the cause of the discrepancy between the calculations 
and experiment. 

PACS numbers: 77.84.-s, 61.50.Ks, 64.70.K-, 61.05.cj 
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I. INTRODUCTION 

Crystals with a general formula ABO3 and the per- 
ovskite structure form a large family whose character- 
istic feature is a tendency to various (ferroelectric, fer- 
roelastic, magnetic, superconducting) phase transitions. 
Barium zirconate BaZrOs, one of the members of this 
family, is a crystal with a high melting point (^2600°C), 
moderate dielectric constant (^40), low dielectric losses, 
and low thermal expansion coefficient. It is widely used 
as a material for microwave ceramic capacitors and for 
inert crucibles used for growth of single crystals of high- 
temperature superconductors. In recent years, the inter- 
est to barium zirconate has revived again because it was 
found that BaZrOs doped with acceptor impurities, in 
particular yttrium, is a high-temperature ionic conduc- 
tor with high proton conductivity, which can be used in 
solid-oxide fuel cells^W Its solid solutions with BaTiOs 
are relaxors with a very high electric-field tuning of the 
dielectric constant;^"— this property can be used in tun- 
able filters, oscillators, phase shifters, etc. It should be 
noted that up to now, the properties of BaZrOa have been 
studied mainly at room and higher temperatures. Only a 
few experiments^ (in addition to the first-principles cal- 
culations^i"— ) were undertaken to study the properties of 
BaZrOs (the ground-state structure, the behavior of the 
dielectric constant, etc.) at low temperatures at which 
most pcrovskite oxides exhibit interesting properties. 

Recent X-ray and neutron diffraction studies^ did not 
found any phase transitions in BaZrOs at T < 300 K, 
so that barium zirconate retains the cubic Pm3m struc- 
ture down to very low temperatures (2 K). These struc- 
tural data, however, disagree with the results of the 
first-principles phonon spectra calculations for BaZrOa 
(Refs . I8l- [ig|l2ll3r ) . in which the appearance of the unsta- 
ble i?25 niode indicates an instability of the cubic struc- 
ture with respect to the rotation of the oxygen octahe- 
dra. In Ref. lid , the equilibrium low-symmetry PI struc- 



ture was found, which results from the doubling the unit 
cell in all directions and the rotation of the ZrOg octa- 
hedra according to Glazer tilt system a~b~c~. Taking 
into account of these rotations was shown to improve 
the agreement between the calculated values of the low- 
temperature dielectric constant with experiment.^'' An 
indirect indication of the possible phase transformation 
in BaZrOs crystals can be its very low thermal expan- 
sion coefficient at T < 300 K4 In order to explain the 
disagreement between the theory and experiment, it was 
suggested that, for some reason, the long-range order in 
the octahedra rotations in these crystals cannot be es- 
tablished even at low temperatures. In particular, it was 
supposed that it may be caused by zero-point lattice vi- 
brationsi^iii 

The aim of this work is to calculate the phonon spec- 
trum and the ground-state structure for BaZrOa from 
first principles and to find experimental evidence for the 
structural instability in this material. We started from 
the assumption that, even if the long-range order in bar- 
ium zirconate is absent for some reason, the structural in- 
stability, if it exists, should be manifested in the change 
of the local structure. To search for these changes, the 
extended X-ray absorption fine structure (EXAFS) stud- 
ies are used. EXAFS is one of the modern powerful 
experimental techniques that enables to obtain detailed 
information on the local structure of crystalsJ^ We ex- 
pected to find the evidence for the structural instability 
in BaZrOa from anomalous behavior of the Debye-Waller 
factors which characterize the magnitude of the local dis- 
placements and the amplitude of thermal vibrations. 



II. CALCULATION AND EXPERIMENTAL 
TECHNIQUES 

In this work, the calculations were performed from first 
principles using the density functional theory. The pseu- 



dopotentials for Ba, Zr, and O atoms used in the cal- 
culation were borrowed from Refs. [isllm The cut-off 
energy of plane waves was 30 Ha (816 eV). The integra- 
tion over the Brillouin zone was performed on the 8x8x8 
Monkhorst-Pack mesh. The lattice parameters and the 
equilibrium atomic positions were determined from the 
condition of decreasing of the residual forces acting on 
the atoms below 5 • 10"'^ Ha/Bohr (0.25 meV/A) while 
the accuracy of the total energy calculation was better 
than 10~^° Ha. To calculate of the phonon spectra, the 
formulas derived from the density functional perturba- 
tion theory and the interpolation techniques, were used. 
The calculated lattice parameter for the cubic BaZrOa 
was ao — 4.1659 A and agreed well with the experimen- 
tal value obtained at 10 K (4.191 A, Ref.l). The slight 
difference of these values was due to a systematic under- 
estimation of the lattice parameter typical for the local 
density approximation (LDA) used in this work. 

The sample of barium zirconate was prepared by the 
solid phase reaction method. The initial components 
were BaCOs and microcrystalline Zr02 obtained by the 
decomposition of ZrOCl2 • 8H2O at 300°C. The compo- 
nents were dried at 600°C, weighed in proper propor- 
tions, ground in acetone and calcined in air at 1100°C 
for 6 h. The resulting powder was ground again and an- 
nealed at 1500°C for 3 h. The single-phase character of 
the sample was confirmed by X-ray diffraction. 

Four EXAFS spectra were recorded using the simulta- 
neous registration of the transmission signal and X-ray 
fluorescence on the KMC-2 station of the BESSY syn- 
chrotron radiation source (the beam energy of 1.7 GeV, 
the beam current up to 300 mA) at the Ba Lni-edge 
(5.247 keV) at 300 K. The intensity of the incident ra- 
diation was measured by an ionization chamber, the in- 
tensity of radiation transmitted through the sample was 
measured with a silicon photodiode, and the intensity of 
the X-ray fluorescence excited in the sample was mea- 
sured using the RONTEC silicon energy-dispersive de- 
tector. 

The Ba Xm-edge was chosen for measurements because 
the expected structural instability, which is accompanied 
by the rotation of the oxygen octahedra, should manifest 
most strongly in the change of the Ba-0 interatomic dis- 
tances, and not in the Zr-0 interatomic distances. This 
is because the Zr-0 distances change very little when 
the octahedra are rotated. In principle, the rotation of 
the octahedra could be also observed at the Zr iiT-edge 
as the changes in the multiple-scattering contributions to 
the EXAFS spectra as proposed in Ref. [l3. 

The EXAFS spectra processing was carried out in the 
traditional way.-iS The spectra were independently pro- 
ceeded, and the results were then averaged. Since, ac- 
cording to our theoretical predictions, the temperature of 
the expected phase transition is below 300 K, the struc- 
ture of BaZrOa at room temperature is cubic, and so the 
structural instability should be manifested as enhanced 
Debye- Waller factors. 
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FIG. 1: Phonon spectrum of BaZrOa in the cubic Pm3m 
phase (a) and the ground-state lA/mcm phase (6). 



TABLE I: The energies of different distorted phases of 
BaZrOa. The energy of the cubic Pm3m phase was taken 
as the energy reference. 



Unstable 
mode 



Glazer 
rotations 



Space group 



Energy 
(meV) 



a a a 


R3c 


-9.17 


a°b-b- 


Imma 


-9.47 


a°a°c- 


lA/mcm 


-10.01 



III. RESULTS OF CALCULATIONS 

The calculated phonon spectrum of BaZrOa with the 
cubic perovskite structure (space group Pm3m) is shown 
in Fig. [1^. It is seen that an unstable phonon with R25 
symmetry (imaginary phonon energies are represented in 
the figure as negative values) appears in this spectrum. 
The eigenvector of the unstable phonon indicates that the 
cubic structure is unstable with respect to the octahedra 
rotations. The results of our calculations are consistent 
with the results of previous calculations,— ~^°'^^i^^ but it 
should be noted that an additional instability at the M 
point, which was observed in Ref. Il2l was not found in 
our phonon spectrum. 

To determine the structure of the ground state, the en- 
ergies of different distorted phases resulting from the con- 
densation of triply degenerate unstable R25 mode were 
calculated. We write a Landau-type expansion of the to- 



tal energy of a crystal in a power series of the amplitudes 
of the distortions, 






(1) 



where -Etot(O) is the total energy of the parent cubic 



phase, and 



tJy, 



and 4>z are the angles of rotation 



around three fourfold axes of the cubic structure. From 
this expansion it can be shown that to find the ground- 
state structure, it is sufficient to calculate the energies of 
phases described by the order parameters (0,0,0), ((/>,0,O), 
and (0,0,0), which correspond to the minimums of -Etot 
for different combinations of signs of the coefficients c 
and e. The specified order parameters correspond to the 
Glazer tilt systems aPaPc', a'^b^b^, and a^a^a", and 
lead to the I A/ mem, Im,ma, and R3c space groups, re- 
spectively. As follows from Table HI the IA/m,cm phase 
has the lowest energy among these phases. This phase 
can be obtained from the cubic PmZm structure by the 
out-of-phase rotation of the oxygen octahedra around one 
of the fourfold axes. The energy of this phase is '^lO meV 
lower than the energy of the cubic phase; this gives an 
estimated phase transition temperature of ^120 K. The 
calculations of the phonon spectrum (Fig. [T})) and elastic 
properties for the lA/mcm phase prove that this phase 
is the ground-state structure. 

It should be noted that among the obtained phases 
there is no PI solution, which was considered as the 
ground-state structure in Ref. 10. This phase corre- 
sponds to the Glazer tilt system a~b~c~, in which the 
rotation angles around three fourfold axes of the cu- 
bic structure are different. The testing of this solution 
showed that the PI structure very slowly relaxes to the 
lA/mcm solution found above (the convergence to this 
result required over 400 iterations). This demonstrates 
that the use of the information on the symmetry of the 
total energy can significantly (by more than 10 times) 
reduce the computation time needed to find the ground- 
state structure and guarantees a correct result. 

The lattice parameters and atomic coordinates in the 
ground-state structure of BaZrOs are given in Table HIl 



IV. EXPERIMENTAL RESULTS 

A typical EXAFS spectrum for the BaZrOa sample ob- 
tained at the Ba Lm-edge at 300 K is shown in Fig. [21 
The data analysis shows that the spectra are consistent 
with a model in which the structure of the sample is 
cubic. The structural parameters (interatomic distances 
and the Debye- Waller factors) for the first three shells are 
given in Table IIIIl It is seen that the structural parame- 
ters obtained using two different measurement methods 
(transmission and fluorescence) are practically the same. 

In Table Hill it should be noted that the Debye- Waller 
factor for the first shell is much higher than the values 




FIG. 2: (Color online) A typical EXAFS spectrum of BaZrOa 
obtained at the Ba Lm-edge at 300 K (thin, black line) and 
the results of its fitting (thick, blue line). 



for the second and third shells (usually the values of this 
factor increase monotonically with increasing interatomic 
distance). 

The values of the Debye- Waller factors are determined 
by two components: (1) static distortions of the struc- 
ture and (2) thermal vibrations. As our measurements 
are performed at a temperature above the temperature 
of the expected phase transition, the static distortions 
of the structure are absent. Furthermore, the optical 
modes responsible for the deformation of the oxygen oc- 
tahedron are high-energy modes in the phonon spectrum 
of BaZrOs and are not excited at 300 K. Therefore, we 
can consider the octahedron ZrOg as a rigid unit. The 
thermal motion of the center of the octahedron relative to 
the Ba atom is characterized by the Debye- Waller factor 
(t|. Since (as result of the rigidity of the octahedron) the 
oxygen atoms are strongly bound to the Zr atoms, the 
only way to explain why af 3> cr| is to take into account 
the rotation of the octahedra. Assuming that in the ab- 
sence of the octahedra rotation the Debye- Waller factors 
for the Ba-0 and Ba-Zr bonds are close (0.008 A^), the 
excess part of the Debye- Waller factor in the first shell. 



Aaf 



ol 



cri, can be associated with the thermal ro- 



tations. Then, using the formula 



Act? = 



,2/)2 



12 



which describes the dependence of Acr^ on the octahe- 
dra rotation angle 6 about the fourfold axis of the cubic 
structure, we can estimate the amplitude of the thermal 
rotations. For the lattice parameter of ag ~ 4.2 A, the 
mean rotation angle is Vff^ ~ 4°. Taking into account 
the results of Ref. |l^ in which the experimental value of 
the rotation angle in SrTiOa at 50 K (well below the tem- 
perature of the structural phase transition) is 2.01°, we 
can conclude the existence of the structural instability in 
BaZrOs. 



TABLE II: Calculated lattice parameters and atomic coordinates in the ground-state structure of BaZrOa 



Phase 



Lattice 
parameters (A) 



Atom 



Wyckoff 
position 



M/mcrn 



a = 5.8722 


Ba 


46 


0.00000 


0.50000 


0.25000 


c = 8.3577 


Zr 


4c 


0.00000 


0.00000 


0.00000 







4a 


0.00000 


0.00000 


0.25000 







8h 


0.22315 


0.72315 


0.00000 



TABLE III: Structural parameters for the first three shells in cubic BaZrOa at 300 K obtained from the EXAFS data. 



Measurement method 


i?i(A) 


al (A=) 


R2 (A) 


ai (A^) 


i?3 (A) 


-1 (A^) 


Fluorescence 
Transmission 


2.916 
2.914 


0.0155 
0.0145 


3.640 
3.635 


0.0079 
0.0068 


4.244 
4.244 


0.0101 
0.0087 


V. 


DISCUSSION 













Now, we discuss why the long-range order in the oc- 
tahedra rotations is absent in BaZrOs at low tempera- 
tures. As was mentioned in Sec. HI one of the reasons 
for this behavior was the suppression of the long-range 
order by zero-point vibrations. In Ref. l2Ci a criterion 
was proposed, which enables, using a set of parameters 
obtained entirely from first principles, to estimate the 
stability of a structure with respect to zero-point vibra- 
tions. According to this criterion, zero-point vibrations 
prevent the system to be localized in one of the local 
minima of the potential if hi'/Eo > 2.419, where v is the 
unstable phonon frequency and Eq is the depth of the 
local minimum of the potential. In our calculations, for 
the unstable R25 mode hi/ = 81. 2i cm~^ « 10 meV and 
£"0 « 10 meV (Table HI); therefore, zero-point vibrations 
should not suppress the structural distortions in barium 
zirconate. The phonon frequency of 79i cm~^ and the 
local minimum depth of 8 meV obtained for BaZrOa in 
Ref. [13 give the ratio hv/EQ ss 1.2, which is close to our 
value. The parameters hv w 27i cm~^ and _Eo = 1.5 meV 
obtained in Ref. [SJ are much lower and give hv/Eo w 2.2, 
which is anyway lower then the critical value. Therefore, 
the explanation proposed in Ref. 8.9, ll| that the sup- 
pression of the long-range order in BaZrOs is caused by 
zero-point vibrations seems unlikely. The fact that the 
influence of zero-point vibrations on the phase transition 
in BaZrOa is weaker than in SrTiOs is not surprising be- 
cause the masses of both metal atoms in barium zirconate 
are larger. 

The theoretical calculations performed in this work en- 
able to propose another explanation for the absence of 
the long-range order in BaZrOs crystals in which the 
structural instability is expected. The closeness of the 
total energies of the R3c, Imma, and M/mcm phases 
(Table ^ , which can result from the condensation of the 
i?25 mode, but do not satisfy the group-subgroup re- 
lation, enables to suggest the appearance of the struc- 



TABLE IV: The minimum and maximum energies of the 
phases resulting from the condensation of rotational R25 and 
Ms modes in some perovskite crystals and the difference of 
these energies. 



Compound 


Unstable 


Space 


Energy 


AE 




mode 


group 


(meV) 


(meV) 


BaZrOs 


R25 


R3c 


-9.17 


0.84 




R25 


I4/mcm 


-10.01 




SrTiOs 


Ms 


P4/mbm 


-9.5 


21.4i^ 




R25 


I4/mcm 


-30.9 




CaTiOa 


Ms 


PA/mhm 


-321 


176i^ 




R25 + Ms 


Pbnm 


-497 




CdTiOg 


R25 


I4/mcm 


-912 


371.2fl 




R25 + Ms 


Pbnm 


-1283 




SrZrOs 


R25 


I4/mcm 


-298 


69 




R25 + Ms 


Pbnm 


-367 





tural glass state when cooling the crystals. The inter- 
nal strains and defects cause the appearance of differ- 
ent phases with close energies at different points of the 
sample, obviously resulting in the absence of the long- 
range order in the octahedra rotations. Moreover, the 
oxygen vacancies, which can be generated in BaZrOa 
quite easily, disturb the three-dimensional connectivity 
of the perovskite structure and can also contribute to 
the structural glass formation. The uniqueness of barium 
zirconate among other crystals with the perovskite struc- 
ture results from the fact that the energy difference be- 
tween the phases with different rotation patterns is only 
0.84 meV (Table HI, which is much less than in other in- 
vestigated crystals (Table irvj) . In our opinion, this small 
energy difference may explain the tendency of BaZrOs to 
the formation of the structural glass. 

Further experiments on BaZrOa, such as studies of 
the temperature dependence of the diffraction line width, 
diffuse scattering, elastic properties, and other physical 



phenomena sensitive to the rotations of the octahedra, 
could provide additional information on the cause of the 
discrepancy between theory and experiment and to con- 
firm or refute our explanation. We think that the low- 
temperature Raman studies of barium zirconate can also 
be very useful; our calculations show that five first-order 
Raman lines located at 16-36, 104-111, 121-128, 350- 
361, and 547-548 cm~^ appear in three possible distorted 
phases (i?3c, Imma, and lA/mcm) with the octahedra 
rotations. 



factor that prevents the establishing of the long-range or- 
der in the octahedra rotations. Very close energies of the 
phases with the different rotation patterns, which can 
result from the condensation of the R25 phonon, sug- 
gest a possible transition to the structural glass state 
when cooling BaZrOs samples. This explains the ap- 
parent contradiction between the theoretical predictions 
and the absence of the long-range order in the octahedra 
rotations observed in the experiment. 



VI. CONCLUSION 
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